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ABSTRACT 

Conventional CCD detectors have two major disadvantages: they are slow to read 
out and they suffer from read noise. These problems combine to make high-speed 
spectroscopy of faint targets the most demanding of astronomical observations. It is 
possible to overcome these weaknesses by using electron-multiplying CCDs (EMC- 
CDs). EMCCDs are conventional frame-transfer CCDs, but with an extended serial 
register containing high- voltage electrodes. An avalanche of secondary electrons is pro- 
duced as the photon-generated electrons are clocked through this register, resulting 
in signal amplification that renders the read noise negligible. Using a combination of 
laboratory measurements with the QUCAM2 EMCCD camera and Monte Carlo mod- 
elling, we show that it is possible to significantly increase the signal-to-noise ratio of an 
observation by using an EMCCD, but only if it is optimised and utilised correctly. We 
also show that even greater gains are possible through the use of photon counting. We 
present a recipe for astronomers to follow when setting up a typical EMCCD obser- 
vation which ensures that maximum signal-to-noise ratio is obtained. We also discuss 
the benefits that EMCCDs would bring if used with the next generation of extremely 
large telescopes. Although we mainly consider the spectroscopic use of EMCCDs, our 
conclusions are equally applicable to imaging. 
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1 INTRODUCTION 

In 2001 a new type of detector was announced: the elec- 
tron multi plying CCD or EM CCD. This was first de- 
scribed by Jerram et al. (20011 ) at E2V Technologies and 
iHvnecek (^OOll ) at Texas Instruments. EMCCDs incor- 
porate an avalanche gain mechanism that renders the 
electronic noise in their readout amplifiers (known as 
read noise) negligible and permits the detection of sin- 
gle photon-generated electrons (or photo -electrons). Whilst 
photon counting in the optical has been possible for 
some time with image tube detectors, suc h as the 
IPCS l|Boksenberg fc Burgess 19721 . 1 Jenkins 19871 ). and with 
avalanche photodiod e-based instruments, such as Optima 
l|Kanbach et al. 20081 ). it has never been available with the 
high quantum efficiency (QE), large format and convenience 
of use of a CCD. 

EMCCDs have generated a lo t of interest in the high 
spatial- resolution community (e.g. iTubbs et al. 2 002). but 
have received much less attention for other astronomical ap- 
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plications. In this paper, we explore how EMCCDs can be 
best exploited for spectroscopy, which is arguably the most 
important and fundamental tool of astronomical research. 
Spectroscopy provides much more information than photom- 
etry, such as the detailed kinematics, chemical abundances 
and physical conditions of astronomical sources. However, 
compared to photometry, where the light from a star or dis- 
tant galaxy is concentrated onto a small region of the de- 
tector, spectroscopy spreads the light from the source across 
the entire length of the detector. The amount of light falling 
onto each pixel of the detector is therefore much lower in 
spectroscopy than in photometry, which makes the reduc- 
tion of read noise much more important and implies that 
EMCCDs should be ideally suited to this application. 



Recognising the potential advantages of EMCCDs for 
astronomical spectroscopy, two separate teams have re- 
cently constructed cameras for this purpose and operated 
them at major observatories: QUCAM2 on the ISIS spec- 
trograph of the 4.2 -m William Herschel Telescope (WHT; 
iTulloch et al. 20091 ) and ULTRASPEC on the EFOSC2 
spectrograph of the 3.5-m New Technology Telescope 



2 S. M. Tulloch and V. S. Dhillon 



jives et al. 20081 . iDhillon et al. 20081 ) 0. These detectors are 
significantly harder to optimise and operate than conven- 
tional CCDs due to the use of higher frame rates and the 
greater visibility of subtle noise sources that can be over- 
looked in a conventional CCD. In fact, incorrect operation 
can actually result in a worse signal-to-noise ratio (SNR) 
due to the presence of multiplication noise in the EMCCD 
(see Section \2. 2 \ . To date, there has been no in-depth anal- 
ysis of the performance of EMCCDs for astronomical spec- 
troscopy presented in the refereed astronomical literature, 
and no guidance given to astronomers on how best to set up 
an EMCCD observation in order to obtain maximum SNR. 
In Section 2, we give a brief review of key EMCCD concepts. 
Section 3 describes the three observing regimes that must be 
considered when using an EMCCD: conventional mode, lin- 
ear mode and photon-counting (PC) mode. The latter mode 
does not lend itself to an analytic treatment, so in Section 
4 we present the results of Monte-Carlo modelling of the 
photon-counting performance of EMCCDs. In Section 5, we 
provide a recipe for obtaining maximum SNR with an EM- 
CCD camera. Finally, in Section 6 we look to the future and 
demonstrate the advantages offered by EMCCDs when used 
for spectroscopy on the proposed 42-m European Extremely 
Large Telescope (E-ELT). 



2 KEY EMCCD CONCEPTS 

Although the principle of operation of an EMCCD is very 
similar to that of a conventional CCD, there are some addi- 
tional features that need to be considered. 



2.1 EMCCD structure 

The structure of an EMCCD (Fig ure 1) has alread y 
been describe d in some de pth by iMackav et al. (200ll ). 
iTulloch (2004h . lMarsh (20081 ) and lives et al. (20081 ). Photo- 
electrons are transferred into a conventional CCD serial reg- 
ister, but before reaching the output amplifier they pass 
through an additional multi-stage register (known as the 
electron-multiplication or EM register) where a high- voltage 
(HV) clock of > 40V produces a multiplication of the photo- 
electrons through a process known as impact ionisation - see 
Figure 2. The EM output amplifier is similar to that found 
in a conventional CCD but is generally faster and hence suf- 
fers from increased read noise. Nevertheless a single photo- 
electron entering the EM register will be amplified to such 
an extent that the read noise is rendered insignificant and 
single photons become clearly visible. Most EMCCDs also 
contain a conventional low-noise secondary amplifier at the 
opposite end of the serial register; use of this output trans- 
forms the EMCCD into a normal CCD. 

Most EMCCDs are of frame-transfer design - see Fig- 
ure 1. Here, half the chip is covered with an opaque light 
shield that defines a storage area. The charge in this storage 
area can be transferred independently of that in the image 
area. This allows an image in the storage area to be read out 
concurrently with the integration of the next image, with 
just a few tens of millisecond dead time between exposures. 
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Figure 1. Schematic structure of an EMCCD, the E2V 
CCD201. Photo-electrons produced in the image area are ver- 
tically clocked downwards, first into the storage area, and then 
into the 1056-pixel serial register. For EMCCD operation, the 
charge is then horizontally clocked leftwards, through the 468- 
pixel extended serial register and into the 604-pixel EM reg- 
ister, before being measured and digitised at the EM output. 
For conventional CCD operation, the charge in the serial reg- 
ister is horizontally clocked rightwards to the normal output. 
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Figure 2. The geometry of the serial and EM registers, show- 
ing how multiplication occurs. The top part of the diagram 
shows a cross-section through the EMCCD structure with the 
electrode phases lying at the surface. Below this are three 
snapshots showing the potential wells and the charge packets 
they contain at key moments (ti,t2,t3) in the clocking process. 
At t3 the photo-electrons undergo avalanche multiplication as 
they fall into the potential well below the HV clock phase. 
Note that this diagram does not show a complete pixel cycle. 



Incorporating frame-transfer architecture into any CCD will 
greatly improve observing efficiency in high frame-rate appli- 
cations where the rea dout time is compar able to the required 
temporal resolution |Dhillon et al. 20071 ). In the case of an 
EMCCD the use of frame-transfer architecture is essential 
otherwise the SNR gains will be nullified by dead time. 



2.2 Multiplication noise 

A single photo-electron entering the EM register can give 
rise to a wide range of output signals. This statistical spread 
constitutes an additional no is e source termed mu ltiplica- 
tion noise (|Hollenhorst 19901 ). iBasden et al. (2003) derive 
the following equation describing the probability p(x) of an 
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Figure 3. Output of an EM register with gA = 100 
in response to a range of inputs from 1 to 5e~ . The 
y-axis shows the probability density function (PDF) of 
the output signal, i.e. the fraction of pixels lying within 
a histogram bin. 



output x from the EM register in response to an input of n 
(an integer) photo-electrons: 



p(x) = 



exp( 



(1) 



This is evaluated for several values of n and with gA the EM 
gain (see Section 2.3) equal to 100, in Figure 3, which shows 
that for an output signal of 300 e~, the input signal could 
have been either 3 or 4e~ with almost equal probability. The 
overall effect is to double the variance of the signal, which 
is statistically equivalent to halving the QE of the camera 
(see Section 3.2). In the photon-noise dominated regime this 
means that conventional CCDs will actually give a higher 
performance. It is in the read-noise dominated regime that 
EMCCDs come into their own, where their lack of read 
noise more than compensates for the effects of multiplica- 
tion noise. Note also that for signal levels where there is a 
low probability of a pixel containing more than one photo- 
electron it is possible to use a photon-counting analysis of 
the image to remove the effect of multiplication noise (see 
Section 3.3). 
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Figure 4. The histogram, plotted on a log e vertical scale, of 
pixels in a weakly-illuminated EMCCD image. The solid line is a 
least-squares linear fit to the photo-electron events lying between 
the two vertical dotted lines. The slope of this fitted line can be 
used to calculate the system gain gs of the camera in e~/ADU. 



a conventional CCD (there are various metho ds, for exam- 
ple the photon transfer curve, [janesick 200 ll ). To measure 
gs, we weakly (< 0.1e~pix _1 ) illuminate an EMCCD with 
a flat field so as to avoid a significant number of pixels con- 
taining more than a single electron. A histogram of such an 
image, with a vertical log scale, is shown in Figure [4] In this 
histogram, the pixels containing photo-electrons lie along a 
curve that is linear except at low values where the effects 
of read noise become dominant. Figure|4]also shows a least- 
squares straight-line fit to the linear part of the histogr am: 
the gradient of the line is equal to — gs l|Tulloch 2004T ). In 
this particular case the camera had a system gain gs = 0.005 
e~/ADU, i.e. a single photo-electron entering the EM regis- 
ter would produce a mean signal of 200 ADU in the output 
image. When discussing noise levels in an EMCCD it is more 
convenient to express this in units of input-referenced photo- 
electrons (ep e ). So in the above example, if the read noise is 
5 ADU this would be quoted as 5 x 0.005 = 0.025e~ e . 



2.3 Gain 

Astronomers typically refer to the gain (or strictly speaking 
system gain, gs) of a CCD camera as the number of photo- 
electrons represented by 1 analogue-to-digital unit (ADU) 
in the raw image, i.e. it has units of e~/ADU. An EMCCD 
camera has another gain parameter that we need to describe: 
the avalanche multiplication gain gA (hereafter referred to 
as the EM gain), and there is a risk of confusion here with 
gs- EM gain is simply a unitless multiplication factor equal 
to the mean number of electrons that exit the EM register 
in response to a single electron input. It is hence related to 
gs by the relation gA = gso/gs, where gso is the system 
gain (in units of e~/ADU) measured with the EM gain set 
to unity. 

To measure the various EMCCD gain parameters, we 
need to first turn off the EM gain by reducing the HV clock 
amplitude to 20V. At this level the EM register will then 
behave as a conventional serial register, i.e. 1 electron in, 1 
electron out. We can now measure gso, just as we would with 



2.4 Clock-induced charge 

Clock-induced charge (CIC) is an important source of noise 
in EMCCDs. Its contribution needs to be minimised. It 
consists of internally-generated electrons produced by clock 
transitions during the readout process. CIC is visible in 
EMCCD bias frames as a scattering of single electron 
events which at first sight are indistinguishable from photo- 
electrons. It is only when a histogram is made of the image 
that they appear different. 

CIC is dependent on a number of factors. The ampli- 
tude of the clo ck swings is relevant, as is the temperature 
l|janesick 200 ll ). At first sight, one might assume that, since 
CIC is proportional to the total number of clock transitions 
a pixel experiences during the readout process, a pixel lying 
far from the readout amplifier should experience a higher 
level of CIC. One would then expect CIC gradients in both 
the horizontal and vertical axes of the image. This would be 
true if the CCD is entirely cleared of charge prior to each 
readout. However, this is never the case. One must consider 
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that prior to each readout the chip has either been flushed in 
a clear operation or read out in a previous exposure. These 
operations leave a 'history' of CIC events in the CCD pixels 
prior to our subsequent measurement readout. The distri- 
bution of these historical events will be higher the closer we 
get to the output amplifier since the CIC charge residing in 
these pixels will have accumulated through a larger num- 
ber of clock transitions than for pixels more distant from 
the amplifier. When these historical events are added to the 
events created in the most recent measurement readout the 
overall effect is that each pixel of the image will have expe- 
rienced the same number of clock transitions regardless of 
its position, and the resulting CIC distribution will be flat. 

CIC is produced by both vertical and horizontal clocks, 
as well as the clocks within the EM register. Vertical CIC 
can be virtually eliminated through the use of non-inverted 
mode clocking (where the clock phases never fall more than 
about 7 volts below the substrate of the CCD, the exact 
value depending on the CCD type). Serial-clock CIC can 
be reduced by using lower clock amplitudes. CIC generated 
within the EM register (described from here on as CICIR) 
is harder to remove since any changes to the EM clock am- 
plitudes produces large changes in EM gain. A well opti- 
mised EMCCD will have its performance limited only by 
CICIR. Other noise sources such as amplifier read noise, 
dark current, image-area and serial-register CIC should all 
have been reduced to an insignificant level with respect to 
C ICIR. This opt imisation process for QUCAM2 is described 
in lTulloch (201(1 ). As an example of this, the histogram of a 
bias image from the QUCAM2 EMCCD camera is compared 
with the histograms of two other images generated using a 
Monte Carlo model (see Section 4) in Figure [5] The first of 
these models consists of CIC originating prior to the EM 
register, the second consists of CIC originating at random 
positions within the EM register. The latter will on average 
experience less multiplication than the former since it will 
pass through fewer stages of the multiplication register, pro- 
ducing a histogram that shows an excess of low value pixels 
(as shown in Figure 0. Various other models were created 
with mixes of the two noise sources. The best fit was found 
to correspond to 85% CICIR and 15% pre-EM-register CIC, 
demonstrating that QUCAM2 has been well optimised (at 
least as far as its CIC performance is concerned with some 
other parameters such as CTE remaining non-optimal.) and 
is dominated by CICIR. Details on the final performance of 
QUCAM2 and the values of some of its more critical param- 
eters can be found in Table[T] Pre-EM register CIC electrons 
should not be confused with dark current generated during 
the readout. These can easily outnumber CIC events if the 
operational temperature is too high, if the CCD controller 
has recently been powered on or if the CCD has recently 
been saturated to beyond full-well capacity. The recovery 
time for these last two cases is approximately 2 hours, and 
an accurate measurement of CIC should not be attempted 
until after such a period. 

CIC electrons generated within the EM register do not 
contribute as much charge to an image as a CIC electron 
generated prior to the register: CICIR has a fractional charge 
when expressed in units of input-referenced photo-electrons. 
If we assume that CICIR is generated randomly throughout 
the EM register then we can calculate the average charge vc 
that it will contribute to an image pixel. The calculation, 
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Figure 5. Histogram (diamonds) of a QUCAM2 bias image com- 
pared with two models: pure in-register CIC and pure pre-EM 
register CIC. A fit to the data, consisting of 85% CIC generated 
within the EM register and 15% generated prior to the EM reg- 
ister is shown as a dotted curve. 



derived in Appendix [X] shows that: 

B c 

Ms a) ' 



vc 



(2) 



where Be is the mean number of CICIR events experienced 
by a pixel during its transit through the EM register. 

2.5 Non-inverted mode operation 

If the vertical clock phases are held more than about 7V 
below substrate then the surface potential of the silicon un- 
derlying the phases becomes pinned at the substrate voltage 
(so-called Inverted-mode operation or IMO). This has im- 
portant consequences for both dark current and CIC. If dur- 
ing the read-out of the image the clock phases never become 
inverted (so called Non-inverted mode operation or NIMO) 
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then the CIC is greatly reduced (in the case of QUCAM2 it 
fell from 0.2e~ pixel -1 to a level that was unmeasurablc 
even after 50 bias frames had been summed). At higher 
operational temperatures NIMO will cause an approximate 
100-fold increase in dark current: something that will negate 
any gains from lowered CIC. The trade-off between CIC and 
dark current does not, however, hold at lower temperatures. 
For QUCAM2 at an operational temperature of 178K, the 
NIMO dark current was 1.5e - pixel - hour - , a value that 
was found to be constant for exposure times of up to 1000s. 
Switching to IMO reduced this to ~0.2e - pixel - hour -1 : a 
value somewhat difficult to measure since it is ~4 times be- 
low the current delivered by cosmic ray events. This demon- 
strates that operation at cryogenic temperatures permits 
NIMO without a loss of performance from high dark cur- 
rent. It should be noted that at higher temperatures, such 
as those experienced by Peltier-cooled CCDs, the situation 
becomes more complex since dark current is no longer con- 
stant with exposure time. This effect was seen during the 
optimisation of QUCAM2 when operated experimentally at 
193K. Using NIMO, the dark current for 60s exposures was 
measured at 12e - pixel -1 hour _1 whereas for 600s exposures 
it was 40e - pixel -1 hour~ 1 . 



3 MODES OF EMCCD OPERATION 

EMCCDs can be utilised in three separate modes, each of- 
fering optimum SNR in certain observational regimes. In 
this section the equations describing the SNR in these three 
modes are shown. 



3.1 Conventional mode 

The SNR obtained through the conventional low-noise am- 
plifier is given by: 



SNR C = 



M 



y/M + is c +D + K + a 2 N 



(3) 



where M is the mean signal per pixel from the source, on 
is the read noise from the conventional amplifier, vc the 
mean CICIR per pixel, D the dark charge and K the charge 
received from sky photons. 

3.2 Linear mode 

In linear mode, the digitised signal from the EM output is 
interpreted as having a linear relationship with the photo- 
electrons, as is usual for a CCD. The SNR obtained through 
the EM output is then given by: 

M 

SNR «" = , , , „ ■ (4) 

y/2.{M + u c +D + K) + (a EM /g A ) 2 

The factor of 2 in the denominator accounts for the multi- 
plication noise ( see section 2.2 ). Th e derivation o f this factor 
can be found in lMarsh f200sh and lTubbs (20031 ). The read 
noise in the EM amplifier a em will typically be tens of elec- 
trons due to its higher bandwidth but its contribution to the 
denominator is rendered negligible by the use of high EM 
gain, gA. High speed means high read noise so high frame 
rate cameras will need higher gains than the more leisurely 
QUCAM2 (1.6s read-out time in EM mode). 
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Figure 6. The effect of the PC threshold value on the detected 
fraction of the signal, the CIC and the read noise in an EMCCD 
image. In this case the read noise uem = 0.025ep e and gj± = 2000. 



3.3 Photon counting mode 

In PC mode we apply a threshold to the image and inter- 
pret any pixels above it as containing a single photo-electron. 
This leads to coincidence losses at higher signal levels where 
there is a significant probability of a pixel receiving two or 
more photo-electrons, but at weak signal-levels where this 
probability is low, PC o peration offers a means of eliminat- 
ing multiplication noise (|Plakhotnik et al. 20061 ) and obtain- 
ing an SNR very close to that of an ideal detector. When 
photon counting we must aim to maximise the fraction of 
genuine photo-electrons that are detected whilst at the same 
time minimising the number of detected CICIR and read 
noise pixels. Figure [S] shows how this can be done. The graph 
shows us that if we set a photon-counting threshold of, say, 
0.1ep e (a value that was later found to be optimum, see 
Section [4. 5p we will detect 90% of photo-electrons but only 
23% of the CICIR. False counts from the read noise will be 
negligible. 

The SNR of an ideal photon-counting detector, includ- 
ing the effects of coincidence losses, is given by: 



SNR BC = 



M 



(5) 



This equation (derived in Appendix|Bj needs to be modified 
to accurately describe a photon-counting EMCCD since it 
makes no allowance for the complex effects of CICIR and 
choice of threshold level. Since the distribution of CICIR 
and photo-electron events are different, the SNR can vary 
greatly depending on the precise choice of threshold. These 
complexities have been explored in detail by modelling (see 
Section|4]) but, in short, the following SNR relation (derived 
in Appendix [Cj is found to hold, assuming a PC threshold 
of 0.1ep e (close to optimum, see Section 4.5): 



0.9M 



V^V ex P [(0-9(M + D + K)/8) + 0.23ln(g A )v c ] 



(6) 



It should be noted that the maximum possible SNR in a sin- 
gle PC frame is w 0.8 (see Figure 11) and it is then necessary 
to average many frames to arrive at a usable image (i.e. one 
with SNR > 3). The number of frames that would need to 
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Figure 7. The relative SNR (compared to an ideal noise-free de- 
tector of the same QE) of each EMCCD mode (conventional, linear 
and PC) is shown over a wide range of illuminations. The solid 
curves show the performance of a detector identical to QUCAM2 
{i>C=0. 013e — , o7v=3.1e — ), the dotted lines shows the performance 
that would be expected from a more highly-optimised detector 
{u c =0.003e~, <rjv=2.6e-). 



be 'blocked' together in this fashion is given the symbol 5 in 
Equation 6. 

3.4 Optimum choice of mode 

Considered from the point of view of maximising the SNR 
per pixel, the choice of readout mode is quite simple. Fig- 
ure 7 shows the range of per-pixel illuminations over which 
each mode offers the best SNR, based on the equations pre- 
sented earlier in this section. A single pixel is, however, not 
usually the same thing as a single wavelength element in a 
reduced spectrum. Many additional factors affect the choice 
of mode, such as plate scale, seeing and sky background, and 
this is explored in greater depth in Section [5] 



4 MODELLING PHOTON-COUNTING 
PERFORMANCE 

Modelling was required for two reasons. First there was no 
equation available describing the output signal distribution 
of pixels affected by CICIR and second to check that the 
assumptions underlying the derivation of Equation [6] (see 
Appendix C) are valid. 

4.1 What was modelled 

Equation Q] describes the output of an EM register for any 
integer number n of photo-electrons. Using the Poisson dis- 
tribution it is possible to calculate the proportion of pixels 



that contain n = 1,2,3.... photo-electrons as a function of 
the mean illumination M. This result can be combined with 
Equation [T] to yield the output distribution of the EM reg- 
ister for any mean input signal level. No equivalent relation 
describing the output distribution of CICIR could be found 
and it is here that Monte Carlo modelling is required (see 
Section 4.2). 

The final output of the model is a pair of 3D vectors. 
The first of these shows the output pixel value distributions 
(i.e. histograms) for a wide range of signal and CICIR levels. 
The second is the associated cumulative distribution func- 
tion (CDF) derived from the histograms in the first 3D vec- 
tor. The CDF is extremely useful since it indicates the mean 
photon counts per-pixel that we can expect for any combina- 
tion of signal and CICIR for any given PC threshold. These 
3D vectors can be visualised as two cubes of normalised his- 
togram values and CDF values. The a;-axis of the cubes are 
labelled with a logarithmic-spaced range of CIC values and 
the y-ax.es with a logarithmic-spaced range of signal values 
(extending up to a maximum of 2.5e~). The z-axes are la- 
belled with pixel values extending up to a maximum of 10e~ e . 
In the case of the CDF cube, the z-axis units can also be 
interpreted as the PC threshold setting and the data values 
as being the mean per pixel-PC signal at that given thresh- 
old. To illustrate what we mean, we show some sample data 
taken from the first of these 3D vectors in Figure El. 

4.2 Modelling of CICIR 

Since no analytical formula describing the distribution of CI- 
CIR events could be found it was necessary to do a Monte 
Carlo model of the EM register. Binomial statistics describe 
processes whose final outcome depends on a series of deci- 
sions each of which has two possible outcomes. It therefore 
applies to the creation of CIC as a pixel is clocked along 
the EM register. Synthetic images were generated contain- 
ing between 1 and 6 CIC events per pixel. It was not really 
necessary to go any higher than 6 CIC events since the bino- 
mial distribution shows that there is an insignificant prob- 
ability of more than 6 events being generated per pixel for 
mean CIC event levels of up to 0.3 per pixel, well beyond the 
useful operational range of an EMCCD (indeed QUCAM2 
gave ~ 0.08 CIC events per pixel). Histograms of these 6 
images were then calculated to yield a set of output CIC 
distributions for integer input, analogous to the distribu- 
tions for photo-electron events given by Equation [T] It was 
then necessary to combine these histograms using the bino- 
mial distribution formula to yield the output distribution of 
the EM register for any mean CIC event level. 

The model was implemented by simulating the transfer 
of charge through an EM register with 604 elements, one 
pixel at a time. At each pixel transfer a dice was thrown 
for each electron in the pixel to decide if a multiplication 
event occurred. An overall EM gain of <7a=2000 was used. 
Six thousand lines were read out in this way to get a good 
statistical sample of pixel values. At the start of the readout 
of each simulated image row, the EM register was charged 
with a single electron per element. This simply amounted 
to initialising the array representing the EM register with 
each element equal to 1. This was then read out, simulat- 
ing the effect of charge amplification, to yield an image with 
width equal to the length of the register. The resulting im- 
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age was then scrambled (i.e. the pixels were reordered in a 
random fashion) and added to its original self to yield an 
image containing an average of 2 CIC events per pixel. This 
scrambling was necessary since the raw images contained 
pixels with values that were approximately proportional to 
their column coordinate, with the pixels in the higher col- 
umn numbers having higher values. Further scramble-plus- 
addition operations were performed to yield images with 3, 
4, 5 and 6 events per pixel. 



4.3 Modelling of realistic EMCCD images 

The distributions of the CICIR (calculated in Section 4.2) 
and the distribution of the photo-electron events were then 
combined, together with read noise, through the use of in- 
termediate model images. Photo-electron events were first 
generated in the image using a random number generator 
that was weighted by the distribution of the EM register 
output. CICIR events were then added to the image in 
the same manner as for the photoelectrons. Finally, read 
noise of <jem = 0.025e~ e was added to every pixel in the 
image. These model images contained a bias region from 
which photo-electrons were excluded. Histograms of the im- 
age and bias regions were then calculated to yield the dis- 
tributions and their CDFs. These CDFs effectively gave the 
mean photon-counting signal from the image and bias areas 
as a function of t the PC threshold. The SNR could then 
be determined using the following equation (derived in Ap- 
pendix [B| : 



SNR pc (i) 



ln[l - CDFi(t)] + ln[l - CDF B (t)\ 



^[l-CDF,(t)]- 



1 



(7) 



where CDFi and CDFb are the image and bias area CDFs, 
respectively. Since the CDFs were calculated over a wide 
range of threshold values it was possible to find the optimum 
threshold value or alternatively just calculate the SNR for 
any given threshold. 



4.4 Testing of the model 

The model output was tested against a stack of 45 QUCAM2 
bias frames of known CIC level, EM gain and system gain. 
The comparison between model and data is shown in Fig- 
ure 8. The agreement is good, although QUCAM2 shows a 
slight excess of low value events compared to that predicted. 
This can be explained by the imperfect charge transfer in the 
EM register which boosts the relative number of low-value 
events, an effect that was not included in the model. 



4.5 Optimum PC threshold 

It is the read noise that sets the lower limit on the photon- 
counting threshold. Gaussian statistics predict that a thresh- 
old set ~ 3<r above this noise gives 1 false count per 1000 
pixels, falling to 1 pixel in 32000 if w e choose a thresh old of 
4cr. Other groups (|Daigle et al. 20061 . lives et al. 20081 ') have 
chosen quite high thresholds (5-5. 5a). This is a good choice 
as long as it is combined with a high EM gain, high enough 
to ensure that the mean level of a photo-electron is at least 
10 times the threshold value. This high ratio between mean 
photo-electron level and read noise permits a threshold to 
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Figure 8. The plot compares the histogram of pixels in a 
genuine QUCAM2 EMCCD image (diamonds) with that of a 
model image (solid line). 



be set low enough to include a majority of photo-electrons. 
There is a limit to how high the EM gain should be pushed, 
though, since it can risk damage to the chip if gain is applied 
during overexposure for long periods. 

Figure 9 shows that the PC threshold needs to be tuned 
depending on the signal level, so that the maximum number 
of genuine photo-electrons are counted and the maximum 
number of CICIR are rejected. Figure 9 also shows that in 
the case of a detector with read noise a em = 0.025e~ e , the 
optimum threshold falls as low as 3.2a em- From a data- 
reduction point of view, using a variable threshold adds 
complexity but may be necessary to extract maximum SNR, 
particularly at low signal levels. One example of this would 
be the measurement of a faint emission line, the peak of 
which would be placed at an optimum signal level through 
a suitable choice of frame rate. Here, the threshold would 
be set low, in the region of 0.1e~ e according to Figure 9. 
The wings of this same line, which may be an order of 
magnitude fainter would then benefit from an increased 
threshold, say ~ 0.25e~ e . The use of an adaptive thresh- 
old would create many side effects, such as noise artifacts 
(the amount of background signal from CIC, sky and dark 
current would be modified depending on threshold setting) , 
so would require extra data reduction effort. Note also that 
non-Gaussian pattern noise is a particular problem in high- 
speed detectors in an observatory environment which may 
require the threshold to be pushed higher than would oth- 
erwise be optimum. The models described so far have as- 
sumed a read noise of 0.025e~ e (equal to that of QUCAM2). 
A simulation was performed of the effect of higher (0.05ep~ e ) 
and lower (0.012ep e ) read noise on the SNR performance 
of an EMCCD. Whilst the higher noise definitely impinges 
on the photon-counting performance by forcing the thresh- 
old higher and giving a lower detected fraction of photo- 
electrons, the lower noise gives very little additional bene- 
fits. This can be explained by the fact that the bulk of the 
CICIR has a distribution that falls between 0-0.05ep e and it 
will dominate any read noise lying within the same range. 

In the specific case of QUCAM2 there is an additional 
reason (apart from staying above the read noise) why the 
threshold must be kept slightly elevated. This is the effect 
of poor charge transfer efficiency (CTE) in the EM register. 
Figure 10 shows the autocorrelation of a low-level flat-field 
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Figure 9. The PC threshold that gives maximum SNR is 
plotted as a function of signal level. Model images were used 
with parameters close to those of QUCAM2. The read noise 
<?EM was 0.025ep e and the multiplication gain gA was 2000. 
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Figure 10. The autocorrelation of a weakly illumi- 
nated QUCAM2 image showing the elongation of the 
single electron events through CTE degradation in the 
EM register. 



which demonstrates the problem. The autocorrelation was 
performed along an axis parallel to the serial register. The 
slight broadening of the autocorrelation peak is indicative of 
less-than-perfect CTE. Using a threshold much below 0.1ep e 
would cause complex effects from multiple counting of single 
electron events due to the slight tail on each event being 

above the thre shold. 

Note that iBasden et al. (2003) have modelled a mul- 
tiple threshold technique that can extend photon-counting 
operation well into the coincidence-loss dominated signal 
regime (see Section 3.3) whilst maintaining high perfor- 
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Figure 11. The SNR of two hypothetical EMCCDs in 
photon-counting mode, one with CICIR Vq equal to that of 
QUCAM2, another with vq equal to that which might be ob- 
tained in a more optimised detector. gA in both cases is 2000 
and the read noise is 0.025ep e . The crosses and stars show the 
predictions of the Monte Carlo model (Equation 7), the solid 
lines show the approximation (Equation [§). The SNR of an 
ideal detector (one where the SNR= y'Signal) is plotted as a 
dot-dash line for comparison. 



is, of course, an approximation and in reality there is a com- 
plex interplay between CIC events and photo-electron events 
in the analogue domain. For example, a photo-electron event 
could be 'helped over' the PC threshold by an accompanying 
CIC event, or a CIC event could be lost by occurring within 
an illuminated pixel. A second simplification would be to 
assume that, as the signal or CIC level increases, the de- 
tected fraction of these events does not change. It is thought 
reasonable to make these approximations since it will only 
seriously fail in the high signal regime where there is a high 
probability of coincidence losses and consequently low SNR 
compared to an ideal detector. The resulting simplified SNR 
equation (Equation [6]) is derived in Appendix [C] This ap- 
proximation was tested against the earlier more comprehen- 
sive model (i.e. that which used Equation 0) for a whole 
range of signals and at two CIC levels. The comparison is 
shown in Figure 11. As can be seen, there is excellent agree- 
ment, justifying our simplifications. 



4.6 Simplification of SNR equation in PC mode 

The method used to calculate SNR in PC mode (Equation]?]) 
is rather complex since it requires the analysis of both im- 
age and bias areas in a large model image. A simpler SNR 
formula that can be applied more generally was therefore 
sought. One simplification would be to consider CIC, read 
noise and photo-electrons separately, i.e. assume that they 
only interact in the digital domain after thresholding. This 



4.7 SNR predictions from model 

The model was used to evaluate the photon-counting SNR 
(SNRp C ) over a range of CIC and signal levels. The results 
are shown in Figure 12, expressed both as a fraction of the 
SNR of an ideal detector SNR ideai and the SNR of an EM- 
CCD operated in linear mode SNRku (as described in Equa- 
tion]?]). The ideal detector is assumed to have the same QE 
as the EMCCD but does not suffer from read noise, coin- 
cidence or threshold losses. The figures reveal the presence 
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Figure 12. Relative SNR contours of a photon-counting EMCCD compared in the left panel to an ideal detector and in the 
right panel to a linear-mode EMCCD. It has been assumed that all the CIC is generated in the EM register, the read noise 
f__M = 0.025e^" e , the threshold is fixed at 0.1ep e and the EM gain, 3^=2000. 



of a photon counting 'sweet-spot', where an SNR in excess 
of 90% of ideal is possible if the CIC can be sufficiently re- 
duced (to around vc = 0.002 e~ pix -1 ). The sweet-spot is 
quite narrow and extends between signals of « 0.07 and 0.2 
e - pix -1 . The plots also show that for signals of less than 
1.2e - pix -1 , photon counting is superior to linear-mode op- 
eration, and this is fairly independent of CIC level. 

Figure 12 shows that the sweet-spot of an EMCCD ac- 
tually covers a very small range of exposure levels, how- 
ever, we can effectively slide the sweet-spot along the per- 
temporal-bin exposure scale to quite high signal levels 
through the use of blocking, i.e. summing together a num- 
ber of frames whose total exposure time equals our required 
temporal resolution. For QUCAM2, if we use the fairly gen- 
erous definition of the sweet-spot as occupying the exposure 
range over which SNR pc > SNR;i„ then this dynamic range 
is about 30:1. This would be like using a normal science 
CCD camera with a 7-bit (16-bit being more usual) ana- 
logue to digital converter (ADC) and could cause problems 
if the spectrum we wish to observe has a set of line intensities 
that exceeds this range. 



5 RECIPE FOR USING AN EMCCD FOR 
ASTRONOMICAL SPECTROSCOPY 

The QUCAM2 camera at the WHT is used here as the ba- 
sis of an example of how to correctly set up an EMCCD to 
maximise the SNR of a spectroscopic observation. The gen- 
eral principles outlined here, however, apply to any EMCCD 
camera. 

QUCAM2 is a relatively slow camera giving a minimum 
frame time in EM mode of 1.6s. This then dictates the high- 
est temporal resolution that is available. The small size of 
the CCD means that when used on the ISIS spectrograph it 



measures only 3.3 arc-minutes in the spatial direction. Full 
frame readout is then generally needed in order to locate 
a suitable comparison star for slit-loss correction and this 
frame rate will therefore be hard to improve upon through 
the use of windowing. The linear EM mode should be con- 
sidered as the default mode about which the observations 
are planned. The reason for this is that it gives an SNR that 
is a constant fraction of an ideal detector for almost all sig- 
nal levels (see Figure 7). The observer will not go far wrong 
by selecting this mode. It may be possible to coax extra 
SNR (as much as 40%) from the observations by switch- 
ing to one of the other modes but if the observation is not 
prepared carefully the data could prove useless. With lin- 
ear mode the observer is guaranteed a practically noise-free 
detector without the dangers of potential coincidence losses 
and worse SNR than with a conventional CCD. 

The first stage in planning the EMCCD observation is 
to refer to Figure 13. This allows us to calculate m, the flux 
per pixel step in wavelength that we can expect from the 
object at our chosen spectral resolution. Using this datum 
we then need to refer to Figure 14 to see how many seconds 
of observation (Tsnri) would be required on our object, us- 
ing an EM detector in linear mode, to reach an SNR=1 in 
the final extracted spectrum. Figure 14 shows the calculated 
times for both dark and bright-sky conditions (new and full 
moon) for both arms of ISIS. Seeing of 0.7" and a slit width 
of 1" has been assumed. Given that the spatial plate-scale 
of QUCAM2 on ISIS is 0.2" pix" 1 , this implies that each 
wavelength element in the final extracted spectrum contains 
~ 5 pixels- worth of sky (assuming the spectrum is extracted 
across ~ 1.5xFWHM pixels in the spatial direction). Once 
we know the time it takes to reach an SNR=1, it is then 
straightforward to calculate the time needed to reach any 
arbitrary SNR, since with linear mode, the SNR is propor- 
tional to the square-root of the observation time. 
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Figure 13. The flux m, that can be expected from each grating on ISIS when used with QUCAM2 as a 
function of source magnitude. i?-band magnitudes as shown for the ij-gratings and B-band magnitudes 
for the B-gratings. Calculated for airmass=l. 
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Figure 14. Tsnri, the observation time on the WHT required to reach an SNR of 1 in the final extracted spectrum with the ISIS 
blue (left panel) and red (right panel) gratings as a function of source brightness. The detector is QUCAM2 operated in linear 
mode. A slit-width of 1", seeing of 0.7" and a spectral extraction over 5 pixels in the spatial direction arc assumed. Observations 
are in the B-band (left) and ij-band (right). 
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The observer can now either play safe and use linear 
mode or explore the possibility of up to 40% higher per- 
formance from either PC or conventional modes. This will 
depend on the mean per-pixel signal (from all sources in- 
cluding the sky) that we can expect during an exposure of 
duration equal to our required temporal resolution, r. This 
is calculated as follows: 

signal pix" 1 = r. |^sky + ^jj , (8) 

where d is the seeing-induced FWHM of the spectrum along 
the spatial axis of the CCD frame, measured in pixels. This 
dictates what spatial-binning factor we later need to use 
when extracting the spectrum. Note that in conventional 
and linear mode, r is defined by the individual frame time, 
whereas in PC mode we divide t into 8 separate frames that 
are later photon counted and summed (in order to observe 
brighter objects without incurring coincidence losses). Now 
that we have an estimation of the per-pixel signal we can use 
Figure 15 to find which mode will give the optimum SNR on 
a per-pixel basis. If we find we can use photon counting then, 
according to Figure 15, the SNR gain will be ~ 25%, which 
we can translate into less telescope time. The SNR increases 
as the square root of the observation time in PC and linear 
mode, since the detector is dominated by noise sources with 
variances that increase linearly with signal. If instead we 
find that we have enough photo-electrons to permit conven- 
tional mode, the reduction in the observation time is harder 
to estimate since the relatively high read noise gives a non- 
linear relation between the square root of the exposure time 
and the SNR. The saving in telescope time from use of con- 
ventional mode could, however, be as much as 50% (relative 
to linear mode) as the per-pixel signal tends to higher and 
higher values and the read noise becomes insignificant rela- 
tive to the Poissonian noise in the sky and target. 



5.1 Worked example of the recipe 

We wish to observe an R=18.5 eclipsing binary star with 
the red arm of ISIS at the highest possible spectral reso- 
lution under dark-sky conditions. The emission line is 4x 
brighter than the underlying continuum. The star under- 
goes an eclipse that lasts 7 minutes that we wish to resolve 
spectroscopically. What SNR can we achieve? Which mode 
should we use? 

To begin this problem, we use Figure 13 to establish 
the signal that we can expect from the emission line. If we 
choose the R1200R grating (resolution ~ 7000) we will re- 
ceive 0.07x4 = 0.28 photo-electrons per wavelength step per 
second. If we observe the object with a temporal resolution 
of 30s we will be able to easily resolve the eclipse. Referring 
to Figure 14 we can then see that this combination of signal 
and spectral resolution will require ~ 7s of observation to 
give an SNR=1 in the final extracted spectrum if we use 
linear mode. Since we can actually observe for 30s, the SNR 
will be equal to ^/30/7 = 2.1. 

We now turn to the choice of observing mode. Assuming 
0.7" seeing, a slit-width of 1" and that the spectrum will be 
extracted over 5 pixels in the spatial direction, we can cal- 
culate that the peak signal per temporal bin per pixel will 
be 0.28 x 30/5 = 1.7 e~. To this we must add the sky signal 
tabulated in Table [Ell Since we are observing in dark time 



and at high resolution this will be a negligible 0.003e~s _1 . 
Referring to Figure 15 we can then immediately rule out 
conventional mode as the SNR would collapse at such a low 
signal level. The best mode would then be PC with 8 some- 
where between 3 and 10. Interpolating the figure we can es- 
timate that 8 w 7 would be optimum. This is feasible since 
the minimum read-time for 7 PC frames is 11.2s with QU- 
CAM2, well below the required temporal resolution of 30s. 
In conclusion, we would obtain the best SNR by observing 
at a frame time of 4.3s, photon counting the raw images 
and then averaging them into groups of 7 to obtain our re- 
quired time resolution. Tuning the exposure time in this way 
to keep the spectral line on the PC sweet-spot affords us a 
~ 25% SNR improvement over linear-mode operation (see 
Figure 15). 

Note that there is an upper-limit to the useful PC block- 
ing factor 8, set by the read noise of the conventional am- 
plifier. Large blocking factors imply low-temporal resolution 
and there comes a point where it becomes favorable to use 
the conventional mode (see Figure 15). The degree of off- 
chip binning f3 used is critical since this adds noise to con- 
ventional mode but not to photon counting mode. 

In the case of QUCAM2 it can be demonstrated (using 
the SNR equations in Section [3| that the maximum useful 
blocking factor is ~ 2Qo%f3. For larger values of 8 the SNR 
that can be obtained with conventional mode then exceeds 
the maximum obtainable with PC operation. 

5.2 Binning 

CCDs can be binned on-chip in a noiseless fashion in both 
axes. Since a spectrum will always be spread by seeing in 
the spatial direction, some degree of binning is often re- 
quired. Some of this can be done on-chip but it is usual to 
do some of it off-chip (i.e. post-readout) during extraction of 
the spectrum. An EMCCD will suffer much less from off-chip 
binning than a conventional detector, which has an effective 
read noise multiplied by the square-root of the binning fac- 
tor. As the off-chip binning factor f3 increases, the balance 
tips ever more in favour of the EMCCD. 

5.3 Phase folding 

For observations of objects that vary on a regular period, 
such as short-period binary stars, we can also consider ex- 
tending our observations over many orbits and then 'phase 
folding' the data. This is an equivalent form of off-chip bin- 
ning. It increases our signal by the folding factor F and per- 
mits us to observe fainter sources or, alternatively, to observe 
the same source at higher time resolution. If the SNR we 
require per wavelength element in our final extracted spec- 
trum is given by SNR re<3 then the number of phase folds F 
we need to use (i.e. the number of orbits we must observe) 
in linear mode is given by: 

^ TSNRI X SNRreq . . 



6 EMCCDS ON LARGE TELESCOPES 

The performance of an EMCCD on an extremely large tele- 
scope and, in particular, the quantitative advantage it might 
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Figure 15. The relative SNR achievable using each mode as a function of the total (i.e. source+sky) per-pixel signal, ft is the 
off-chip binning factor. <5 is the photon counting blocking factor (the number of thrcsholdcd PC frames that are summed within 
each temporal bin). The read noise and CIC levels experienced with QUCAM2 are assumed. 
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Figure 16. The temporal resolution versus source magnitude for two seeing conditions on an ISIS-typc E-ELT instrument. The 
curves indicate the observation time required to reach an SNR of 1 in each wavelength clement of the final extracted spectrum 
using a normal (non-FT) CCD, a FT CCD, an EMCCD in linear mode and an ideal detector. Off-chip spatial binning of x2 and 
<tjv = 2.6c~ arc assumed. The grating is a VPH equivalent of the ISIS R1200R grating and we observe in the J?-band during dark 
time. The slit width is assumed to be 1.5x the seeing FWHM. The curve corresponding to the non-FT CCD assumes a 2k x 4k 
detector (similar to that currently used on ISIS) with full-frame read out. An ideal detector is one with the same QE as a normal 
CCD but which has zero read noise. 
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give over a conventional CCD camera is explored in this 
section. Signal fluxes are calculated assuming a hypotheti- 
cal ISIS-type instrument on the E-ELT (42-m aperture) us- 
ing higher efficiency volume-phase holographic (VPH) grat- 
ings. These typically offer a 30% increase in throughput 
over conventional diffraction gratings. Sky backgrounds from 
Paranal (approximately 20 km dista nt from the future E- 
ELT site) are assumed (|Patat 20031 '). The QE of current 
EMCCDs are already very high (peaking at > 90%) and 
any future EMCCD is unlikely to be significantly different. 
The temporal resolution as a function of source magnitude is 
calculated assuming that an SNR of 1 is required in each ele- 
ment of the final reduced spectrum. Equations [3] and [4] were 
used for these calculations. The results are shown in Fig- 
ure 16, where linear and conventional modes are compared 
with an ideal detector. Two seeing conditions are assumed: 
natural median seeing (0.6") and that which may be ob- 
tained in the ii-band with adaptive optics (AO) correction 
(0.1"). With natural seeing and when observing very faint 
sources, long integrations are required. Under these condi- 
tions, the sky contribution becomes so high that the use of 
EMCCDs actually degrades the SNR through the effects of 
multiplication noise. If instead AO is used, the sky contri- 
bution is considerably reduced and the EMCCD operated 
in linear mode gives an advantage right across the range of 
source intensities explored in the plot. 

6.1 Need for high frame rates 

Higher frame rates allow us to observe at higher time resolu- 
tion, of course, but in the case of PC operation it also means 
we can avoid coincidence losses when observing brighter 
sources and cope with the higher sky backgrounds we can ex- 
pect from the E-ELT's larger collecting area. We always need 
to ensure that the per-pixel signal remains in the sweet-spot 
(< 0.2e~, see Figure 7), so even in dark time and at high 
spectral dispersion we would need (referring to Table IE2() 
to operate at > 0.5Hz in the blue and >2Hz in the _R-band 
if we wish to photon count on the E-ELT (natural median 
seeing is assumed). 

6.2 Need for frame transfer design 

Frame-transfer design reduces the dead time between expo- 
sures to a few tens of milliseconds: the time it takes to move 
the image into the storage area. For a normal mechanically- 
shuttered non-FT camera, the dead time is equal to the 
read-out time of the CCD. This can have a dramatic effect 
on observing efficiency. For example, assuming a scientific 
camera is used to take blocks of 20 minute exposures for 10 
hours a night, then with a read-out time of 60s this dead time 
could amount to more than two weeks of telescope time over 
the course of a year! In high time-resolution spectroscopy 
with much shorter exposure times, the dead-time losses are 
proportionately more extreme, resulting in unfeasibly low 
observing efficiencies; it is here that an FT CCD can give 
massive gains. The performance of a conventional CCD with 
non-FT architecture and 30s dead time between exposures 
is shown plotted in Figure 16. Any future CCD used on the 
E-ELT, especially if destined for high time-resolution appli- 
cations, be it an EMCCD or otherwise, should incorporate 
FT architecture. 



6.3 An EMCCD for the E-ELT 

Current EMCCDs are rather small (lk x lk pixels). For spec- 
troscopy this is not a good format since it limits our spec- 
tral range and also the availability of comparison stars near 
the target. Any EMCCD designed for use with the E-ELT 
therefore needs to be physically larger. It will require mul- 
tiple outputs to achieve higher frame rates, so as to permit 
high time-resolution and allow photon-counting mode oper- 
ation. A possible geometry for such an EMCCD is shown in 
Figure 17. Here a monolithic 4k x 2k image area is proposed. 
Tapered storage areas are positioned both above and below 
the image area. Their shape provides space for multiple EM 
registers. This is similar to the architecture of the CCD220, 
a smaller format, multiple output EMC CD developed for 
wavefront sensing (jFeautrier et al. 20081 ). In principle, the 
storage areas can be shrunk by reducing the height (in the 
axis perpendicular to the serial register) of their pixels. This 
results in more efficient use of the available silicon. It would 
be at the expense of reduced full-well capacity, but this is 
of little concern in the low-signal regime where EMCCDs 
are best applied. Each storage area can be read out through 
either a conventional amplifier capable of giving 2.6e~ noise 
at 200kpix s _1 or a higher-noise (20-30e~) EM amplifier ca- 
pable of running at lOMpix s _1 . These amplifiers could be 
combined so that conventional mode would be selected sim- 
ply by setting the EM gain to unity, but having two separate 
amplifiers allows more design flexibility with one amplifier 
being optimised for speed and the other optimised for low 
noise. The target CICIR level should be around 0.003e~. 
Non-inverted mode operation and liquid nitrogen cooling 
would ensure that dark current and CIC generated prior to 
the EM register are kept at negligibly low levels. Full-frame 
readout would be possible in 5s through the conventional 
amplifiers and in 0.1s through the EM amplifiers. None of 
these design parameters, taken individually, exceed those of 
current scientific cameras. We already know from discussions 
with E2V that such a device is feasible to manufacture, al- 
though the large pin-count may require the EM and conven- 
tional amplifiers to be combined. The spectral axis would lie 
horizontally in Figure 17. This would maximise the spectral 
coverage and also allow efficient on-chip spatial binning in 
order to reduce the effects of CICIR and read noise, and to 
reduce the read time in close proportion to the spatial bin- 
ning factor. Such a CCD would have a wide application in 
the field of high time-resolution spectroscopy. In Figure 16 
it has already been shown that simply switching to an FT 
design without implementing an EM register can give a huge 
advantage. Considering the economics of the E-ELT, the ex- 
tra expense required to implement FT architecture (which 
approximately doubles the area of silicon required) can be 
easily justified. The incremental expense of then adding an 
EM register (estimated at 10-20% by E2V) will give excel- 
lent value for money given the further savings in telescope 
time that it can provide. 



7 CONCLUSIONS 

We have shown that EMCCDs are almost perfect detectors 
for optical spectroscopy. They have close to 100% QE, virtu- 
ally no read noise, large formats, linear response, negligible 
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from. The biggest gains, however, will be for astronomers 
doing read noise limited spectroscopy, due to the fact that 
the read noise will also be zero. It is our firm belief, there- 
fore, that once these large-format EMCCDs become avail- 
able, they will become the detector of choice on the world's 
major spectrographs including those to be built for the E- 
ELT. 
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CE CE CE CE 

Figure 17. A possible design for a future EMCCD for spec- 
trographic use on the E-ELT. Conventional amplifiers are la- 
belled C, EM amplifiers are labelled E. The shaded regions 
are the storage areas. Each lk pixel x lk pixel block has its 
own conventional and EM amplifiers, permitting faster paral- 
lel readout. This monolithic design is buttablc along two edges 
allowing it to be extended horizontally for increased spectral 
coverage. 



dead time and are relatively inexpensive. Due to multipli- 
cation noise and CIC, however, they are not as straightfor- 
ward to use as conventional CCDs, and care must be taken 
to ensure that they are operated in the correct mode: con- 
ventional, linear or photon counting. If EMCCDs are used 
correctly, it is possible to gain orders-of-magnitude improve- 
ment in SNR compared to conventional CCDs at low light 
levels (see Figure 7) , but using the wrong mode can result in 
an orders-of-magnitude reduction in SNR. With the aid of 
Monte-Carlo modelling of an EMCCD, we provide deriva- 
tions of the SNR equations for each EMCCD mode, and 
present a recipe for astronomers to assist in determining the 
optimum EMCCD mode for a given observation. 

The EMCCD used in QUCAM2 and ULTRASPEC is 
the E2V CCD201-20 detector, which at lkxlk pixels (each 
of 13/imxl3/im) is the largest commercially available EM- 
CCD. However, this is still substantially smaller than the 
conventional 4kx2k CCDs found on major optical spectro- 
graphs. As a consequence, using the CCD201 results in the 
loss of approximately three-quarters of the wavelength cov- 
erage and one half of the spatial coverage provided by typical 
astronomical spectrographs. This factor of eight loss of de- 
tector area is a heavy price to pay, even for the huge SNR 
gains of an EMCCD. We therefore present a concept for a 
large-format EMCCD which can optimally sample the focal 
plane of the world's major spectrographs. It is important 
to emphasize that EMCCDs are identical to conventional 
CCDs in almost every respect (architecture, performance, 
read-out electronics), except for the fact that they have ef- 
fectively zero read noise. This means that astronomers who 
do not have read noise limited observations lose nothing by 
using an EMCCD instead of a conventional CCD; in fact, 
we show that there will be significant efficiency gains (equiv- 
alent to approximately two weeks of time per year per tele- 
scope) because the frame-transfer format of EMCCDs re- 
sults in essentially zero dead time compared to the tens-of- 
seconds dead time that conventional non-FT CCDs suffer 
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APPENDIX A: FRACTIONAL CHARGE OF 
CICIR 

A CIC electron originating within the EM register will ef- 
fectively have a fractional charge whose value is equal to the 
average charge qo of a single photo-electron charge packet 
during its transit through the EM register. The instanta- 
neous charge qo of a pixel within the EM register that orig- 
inated as a single photo-electron at the register input is: 



(1+5 



(Al) 



where x is the position within EM register and p the per- 
transfer multiplication probability. The mean value, qo, of 
this pixel during its EM register transit can then be obtained 
by integrating this function over the length of the register 
and then dividing by the number of stages S. If we refer- 
ence this charge to an equivalent signal at the input to the 
register, such that qi — qo/g A , we get: 

(A2) 



1 f a 

II = -n / (1 +pYdx. 

b 9A J x=1 



This integral by the standard solution 
1 



qi = 



(1+P)"]f =1 , 



(A3) 



1 



(A4) 



Sg A \n(l+p) 
which gives the result 

In g A In g A ' 

Equation IA4I now allows us to calculate the per transfer 
probability of CIC pc in the EM register from a knowledge 
of the gain g A and the mean CIC charge in the bias vc- 
The generation of this charge is a binomial process so Be 
the total number of CIC events in a pixel exiting the EM 
register is given by 



B c = pcS. 



(A5) 



Since each of these CIC events will contain an average charge 
of qi , the mean per pixel CIC charge vc is given by 

vc = qiB c , (A6) 

substituting qi from Equation I A5I and rearranging, we get 

\n(g A )v c 



PC 



S 



(A7) 



We can then show how the mean CIC charge in a pixel vc 
(relevant for linear-mode operation) and the mean number 



of CIC events in a pixel Be (relevant for PC operation) are 
related as follows: 

Be 
vc 



In(g/ 



(A8) 



APPENDIX B: SNR OF AN IDEAL PHOTON 
COUNTER 

We present a full derivation of SNR pc , the SNR in a photon- 
counting detector. This differs from that of an ideal detector 
due to the effects of coincidence losses. To the best of our 
knowledge this derivation has not appeared before in the 
astronomical literature. 

Let N be the mean illumination in photo-electrons 
per pixel and n the mean photon-counted signal per pixel, 
i.e. the fraction of pixels that contain one or more photo- 
electrons. Poisson statistics tells us that 



n = 1 



therefore: 



N = -ln(l-n). 



(Bl) 



(B2) 



The noise in a photon counted frame can be derived straight- 
forwardly by considering that only two pixel values are pos- 
sible: and 1. Pixels containing will have a variance of 
N, those containing 1 will have a variance of TV — 1. Know- 
ing the fraction of pixels containing each of these two values 
then allows us to combine the variances in quadrature to 
yield a pc , the mis noise: 



= y/[e- N N 2 + (1 - e~ N )(N - l) 2 



(B3) 



(B4) 



The photon-counted images must then be processed to re- 
move the effects of coincidence losses. This is done after the 
component frames within each temporal bin have been aver- 
aged to yield a mean value for n for each pixel. The original 
mean signal N prior to coincidence losses is then recovered 
by using Equation IB2I Although coincidence loss tends to 
produce a saturation and a smoothing of the image struc- 
ture, the overall effect is to add a great deal of noise to the 
observation and for this reason we must avoid a photon- 
counting detector entering the coincidence loss regime. The 
amount of extra noise generated can be calculated by con- 
sidering the change dN in TV produced by a small change 
dn in n. The noise in the final coincidence-corrected pixel 
will then be equal to that in the unprocessed average pixel 
multiplied by dN/dn. From Equation IB2I we get 

N + dN = - ln[l - (n + dn)]. (B5) 

This standard differential is then solved to yield 

§ = (l-»)- 1 . (B6) 
an 

Substituting N for n using Equation IB II we get 

^ = e N - (B7) 
dn 

We then multiply the uncorrected noise given in Equa- 
tion IB4I by this factor to yield the noise in the final 
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coincidence-loss-corrected PC image. SNR pc is then given 
by: 

N 



SNR DC = 



(B8) 



s/e N - 1 

This can be expressed in units of n, using Equation IB21 
which is more useful since it is n that we actually measure 
from our images: 



SNR pc = 



■ln(l-n) 



va-^)- 1 



(B9) 



APPENDIX C: SNR OF A 
PHOTON-COUNTING EMCCD 

We have already shown in Appendix [B] that the SNR of an 
ideal photon counter is: 



SNR BC = 



N 



(CI) 



with N representing the signal per pixel. This basic equa- 
tion is now altered to a more realistic form to include the 
noise sources found in an EMCCD. Certain approximations 
are made during this process. The validity of these approx- 
imations have been verified by the Monte Carlo modelling. 

So to begin with, we replace N in the numerator with 
the detected fraction of photo-electrons at our chosen thresh- 
old and we replace N in the denominator with the detected 
signal plus the detected CIC. In Figure 9 we have already 
shown that a threshold of 0.1e pe is close to optimum and 
that at this level 90% of photo-electrons and 23% of CICIR 
will be detected. We then get: 



SNR pc = 



0.9M 



^exp(0.9M + 0.23Bc) 



(C2) 



where M is the signal per temporal bin and Be the number 
of CICIR events per pixel. We now need to consider that any 
photon counting observation will require the blocking of 5 
separate images if we are to achieve a usable SNR (the SNR 
of a single PC image with the exposure level lying within 
the sweet-spot is w 0.4, see Figure 11). We then get: 



SNR D 



0.9M 



VS^exp(0.9M/S + 0.2Wc) 



(C3) 



Next we need to include other noise sources such as sky K 
and dark charge D (units of e~ pix" 1 ). Since these are in- 
distinguishable from photo-electrons they will have the same 
detected fraction for a given PC threshold. We also need to 
express the CICIR in terms of vc (see Equation I A8|) . i.e. the 
mean per-pixel charge that CICIR contributes to the image. 
This is a parameter that can be measured directly from the 
bias frames of an EMCCD camera. So we get: 

0.9M 

+ D + K)/S) + 0.231ii( 5 a)^c] -T 

(C4) 

Note that vq the CICIR is multiplied by a factor of In {gjC} 
in the denominator of Equation lC4l This is a consequence of 
the fractional charge of a CICIR electron (see Appendix [A)) . 

Note also that the read noise has been entirely ignored 
in this simplified description. This is fair since the thresh- 
old was set well above the noise, resulting in very few false 



counts. The read noise has a secondary influence, however, 
since it causes an effective blurring of the threshold level. 
This 'fuzzy-threshold' must add some noise to the images 
since it can make all the difference as to whether an event ly- 
ing close to the threshold is counted or not. The close fit be- 
tween the approximation and the comprehensive model (see 
Figure 11) would, however, indicate that this noise source is 
not significant. 



APPENDIX D: MODEL VARIABLES 

The data cube describing the count-rate from an EMCCD 
in PC-mode as a function of the mean illumination, the 
mean number of CIC events per pixel and the thresh- 
old setting, is available as an IDL .sav file, together with 
descriptions of the variables it contains, at this address: 
http:/ /www. qucam.com/emccd/Histogram. html 
It is included for those wishing to check our results or to 
perform their own experiments. Figure Dl shows a selec- 
tion of histograms extracted from the data cube by way of 
example. 



1 c° r 




Figure Dl. Some histograms generated by the EMCCD 
model, described in Section 4, are plotted here. These are 
taken from the IDL variable outputMIXhist. In the upper panel 
are a set of histograms with the signal level varying from 
0.001 to 2.5e~ pix" 1 and with the CIC fixed at v c = O.OOle- 
pix -1 . The lower plot has the same range of signals but with 
vc = 0.03c~ pix -1 . The read noise was 0.025ep e in both 
cases. 



Astronomical spectroscopy with EMCCDs. 



18 S. M. Tulloch and V. S. Dhillon 



Table El. WHT+ISIS sky backgrounds with QUCAM2. Units: 
photo-electrons pix - 1 s — Slit width= 1". From the ING web 
pages. 



Grating 


Bright 


Dark 


R316R 


0.12 


0.014 


R600R 


0.05 


0.006 


R1200R 


0.023 


0.003 


R300B 


0.07 


0.004 


R600B 


0.03 


0.002 


R1200B 


0.016 


0.001 



Table E2. Estimated sky backgrounds for an E-ELT+ISIS-type 
instrument with VPH gratings. The plate scale and the detector 
QE is assumed to be the same as for QUCAM2 on ISIS. Units: 
photo-electrons pix _1 s _1 . Slit width= 1". 



Grating Bright Dark 



Red 316 


15 


1.6 


Red 600 


6.0 


0.7 


Red 1200 


2.7 


0.4 


Blue 300 


8.2 


0.47 


Blue 600 


3.6 


0.23 


Blue 1200 


1.9 


0.12 
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This paper has been typeset from a TffjX/ WFftK file prepared 
by the author. 



